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Неавтономне виявлення викидпв 1 очищення 
сигналив мон1торингу внутрииньочерепного тиску 


Сигнал мониторинга внутричерепного давления (ВЧД), полученный в отделениях нейрореанимации, 
часто содержит большое количество шумов и выбросов. Эти артефакты не только непосредственно 
приводят к ложным тревогам в автоматических аварийных системах мониторинго-диагностических 
комплексов для управляемой терапии больных, но и сильно загрязняют особенности основного сигнала, 
что делает невозможным точное прогнозирование вторичных повреждений головного мозга, вызванных 
внутричерепной гипертензией. В статье предлагается эффективный онлайновый двухэтапный метод 
очистки физиологических сигналов, основанный на идентификаторе Хампеля и калмановской фильтрации. 
Вначале измеряемые ВЧД сигналы проходят предварительную обработку, в которой осуществляется 
идентификация структуры данных сигнала, оценивается уровень шума, и удаляются выбросы с помощью 
робастного фильтра Хампеля. Точки выбросов заменяются значением медианы в этих точках. Затем 
производится коррекция этих точек, удаление шума и оценка спектра сигнала на основе адаптивной 
авторегрессионной (ААВ) модели с помощью фильтра Калмана и сопряженного с ним КТ$ (КаисВ- 
Тапэ-5еБе!) сглаживающего фильтра. Настраиваемыми параметрами предлагаемого метода фильтрации 
являются полуширина скользящего окна, пороговая величина для обнаружения выбросов и порядок 
ААК модели 

Ключевые слова: фильтр Хампеля, калмановская фильтрация, адаптивная авторегрессионная 
модель, выбросы, внутричерепное давление. 
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Сигнал монторингу внутрипньочерепного тиску (ВЧД), здобутий у вддленнях нейрореанмацй, часто 
мстить велику юмльк!сть шумв 1 викидв. Ц! артефакти не лише безпосередньо викликають помилков1 
тривоги в автоматичних авар1йних системах мониторинго-д1агностичних комплексов для керованот 
терапй хворих, але й сильно забруднюють особливост! основного сигналу, що унеможливлюе точне 
прогнозування вторинного пошкодження головного мозку, викликаного внутрииньочерепною гтертензиею. 
У статт! пропонуеться ефективний двоетапний метод онлайнового очищення ф!з1олойчних сигналв, 
заснований на 1дентиф!катор! Хампеля 1 калмановсько! ф1льтрацй. Спочатку вимрюван! ВЧД сигнали 
проходять попередне оброблення, шд час якого здснюеться 1дентиф!кащя структури даних сигналу, 
ощнюеться р!вень шуму 1 усуваються викиди за допомогою робастного ф!льтра Хампеля. Точки викид!в 
замтнюються значенням мед1ани в цих точках. Пот!м в1дбуваеться корекщя цих даних, усування шуму 
1 ощнка спектра сигналу на основ! адаптивно! авторегрес!йно! (ААВ) модел! за допомогою фультра 
Калмана 1 зв’язаного з ним КТ$ (Ваисв-Типе-561ебе]) згладжувального фультру. Регульованими параметрами 
пропонованого методу ф1льтраци е швширина ковзного в!кна, порогова величина для виявлення викид!в 
1 порядок ААК модел:. 

Ключов! слова: флльтр Хампеля, калмановська фультращя, адаптивна авторегресйна модель, 
викиди, внутрииньочерепний тиск. 


Введение 


Подавляющее число исследований по цифровой обработке данных относится к 
сигналам техногенного происхождения, в которых решаются задачи помехоустойчи- 
вости и фильтрации зашумленных сигналов, формируемых различными средствами 
связи и вещания. В последние десятилетия, особенно за рубежом, стали интенсивно 
развиваться компьютерные методы обработки физиологических сигналов в медицинских 
приложениях. Наиболее развитыми в этом отношении являются методы обработки 
сигналов электрокардиограмм, электроэнцефалограмм, ультразвуковой диагностики. 
На основе современных информационных технологий возникли методики рентге- 
новской, магнито-резонансной, позитронно-эмиссионной томографий. Анализ физи- 
ологических сигналов оказался гораздо сложнее, чем техногенных. Это обусловлено 
двумя обстоятельствами: физиологическим процессам свойственны сигналы с низкими 
уровнями амплитуды и эти сигналы по своей природе являются нестационарными и 
нелинейными. Более того, эти сигналы загрязнены в большой степени стохастическими 
компонентами и случайными выбросами, произведенными другими органами тела, 
мускульными сокращениями и дыханием. 

Выбросами в физиологических временных сигналах являются аномальные данные 
наблюдения, которые существенно отклоняются от большинства данных наблюдения 
и кажутся изолированными во временном ряду. Они могут возникать в результате 
шума измерительных детекторов, сбоя измерения и передачи данных в системах 
управления, внешних возмущений на измерительный процесс, например, из-за наводки 
на соединительные провода и детектор, а также из-за человеческого фактора при 
мониторинге диагностических биомедицинских сигналов у пациентов. Проводить 
компьютерный анализ загрязненных выбросами сигналов, на основе которого ставится 
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диагноз, бесполезно, потому что выбросы могут привести к ошибке конкретизации 
эмпирического моделирования, тенденциозной оценке параметров модели и к непра- 
вильным аналитическим результатам. Слабый физиологический сигнал, например, 
такой, как сигнал внутричерепного давления (ВЧД), главным образом погружен в 
артефакты движения, сопровождающиеся выбросами. Прежде чем делать выбор под- 
ходящего алгоритма цифровой обработки этого сигнала, необходимо использовать 
надлежащий метод предварительной обработки, чтобы ослабить мощность выбросов — 
помехи движения, сохранив при этом данные, искаженные только флуктуационной 
помехой. С этой точки зрения, предварительная обработка данных отличается от 
фильтрования данных. Фильтрование данных изменяет их структуру, не только удаляя 
выбросы, но также и уменьшая вариации данных. Предварительная обработка данных 
является в основном более трудной задачей, чем фильтрование, потому что она одно- 
временно требует идентификации структуры данных и оценки уровня шума, чтобы 
удалить выбросы и сохранить «хорошие» данные. 

Для обнаружения и удаления выбросов из обрабатываемых данных временного 
ряда существует, по-видимому, три базовых фильтра, которые в недавних работах 
получили дальнейшую модернизацию. Это медианный фильтр, близко связанный с 
ним и почти также прост — фильтр Хампеля (Натре! НКег) [1] и фильтр-очиститель, 
разработанный Мартином и Томсоном (Мага апа Твотзоп, МТ ЯКег-сеапет) [2]. 
Если два первых фильтра работают в неавтономном (онлайновом) режиме, то МП 
фильтр-очиститель — по существу автономная схема, которая основана на предвари- 
тельно оцененной модели временного ряда. Этот МП фильтр-очиститель состоит из 
трех шагов. Вначале модель обработки данных временного ряда подгоняется линейной 
авторегресионной АК(р) моделью низкого порядка, используя робастные обобщенные 
М-оценочные функции. Учитывая модель аддитивного выброса (АО), был разработан 
алгоритм фильтра в виде пространства состояний для последующего использования 
фильтра Калмана, в котором, с целью уменьшения выбросов, к последовательности 
обновлений применена Ф-функция. Практически, Ф-функция в основном выбирается 
как правило «трех-сигм» («Зо ед ге»). Наконец, калмановская фильтрация гауссовой 
части данных временного ряда используется для последующей спектральной оценки. 
Конечно, этот МП фильтр-очиститель сложен для практического применения, к тому 
же пороговая точка фильтра низка, когда порядок модели становится высоким. 

В статье [3] авторы предложили модифицированный МП фильтр-очиститель, 
основанный на скользящем окне, которое может захватить динамические изменения 
данных временного ряда в процессе онлайнового режима. У предложенного фильтра- 
очистителя также есть фиксированная пороговая точка для любого порядка модели. 

Эти три базовых фильтра проходили тестирования на идентификацию трех раз- 
личных динамических систем, в одной из которых обрабатывался реальный, зашум- 
ленный выбросами, сигнал команды управления высотой и углом наклона полета 
вертолета и две другие модельные системы - линейная модель и нелинейная модель 
Волтерры второго порядка с аддитивными выбросами и гауссовыми шумами, для 
которых известны точные результаты [4]. Как оказалось, выбросы действительно 
неблагоприятно влияют на результаты анализа обрабатываемых данных или эмпири- 
ческого моделирования, для эффективной очистки фильтры должны, в основном, 
быть нелинейными, хотя линейные фильтры действительно уменьшают выбросы, но 
они не устраняют их и вводят нежелательное уширение. И наконец, по мнению автора [4], 
не существует эффективного универсального фильтра очистки данных. Это наиболее 
наглядно на его примере с медианным фильтром, который был наиболее эффективным 
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для вертолетных данных, но наименее эффективным для модельной нелинейной сис- 
темы Вольтерры. Точно так же, МП фильтр-очиститель оказался самым эффективным 
для линейного примера моделирования и неэффективным в смоделированном примере 
Волтерры, и полностью неподходящим к вертолетным данным, поскольку номинальная 
вариация в вертолетной последовательности данных далека от гауссовой. 

Целью данной работы является построение основанного на идентификаторе 
Хампеля, простого робастного фильтра со скользящим окном для удаления выбросов 
в неавтономном режиме в комбинации с фильтром Калмана и сопряженным с ним 
ВТ$ (Каисв-Типз-5 ее!) сглаживающим фильтром для прогноза и коррекции ори- 
гинальных особенностей основного сигнала внутричерепного давления (ВЧД). 


Методы поиска и замены выбросов 
1 Робастное оценивание выбросов — фильтр Хампеля 


Большинство методов обнаружения выброса основано на базовом предположении, 
что случайные переменные независимы и равномерно распределены, т.е. распределение 
вероятности каждой случайной переменной то же самое, что и у других, и все они 
взаимно независимы. Двумя самыми важными статистическими данными для анализа 
данных в присутствии выбросов являются местоположение (например, среднее) данных 
сигнала и их разброс (дисперсия, ковариация). Классическим методом оценивания 
параметров модели является метод наименьших квадратов, однако на практике не 
всегда выполняется гауссово распределение измеряемых данных и если равномерность 
распределения нарушается даже одним выбросом, его влияние, как на выборочное 


среднее Х, так и на дисперсию 5? выборки Х, = {х,}”, может значительно исказить 


статистические параметры выборки. В связи с этим возникает надобность в построении 
робастных методов оценивания, которые должны обнаружить выбросы и удалить их 
либо заменить на интерполированные значения. 

Для измерения устойчивости М-оценочных функций в зависимости от влияния 
на них выбросов Хампель [5] ввел понятие пороговой точки. В упрощенном виде 
пороговой точкой является наименьшее процентное соотношение загрязнения данных 
выбросами, которое может привести к неограниченному значению ошибки. В целом, 
чем выше пороговая точка, тем более робастная оценочная функция. Для М-оценок 
пороговая точка не должна превышать 1/М№ , где М — размерность параметрического 
пространства [5], поэтому при даже умеренном числе параметров значение пороговой 
точки становится очень малым и в случае оценок наименьших квадратов она стремится 
к нулю. В общем случае робастной оценки выбросами можно считать те величины, 
которые удовлетворяют критерию: 


| —жо) >16, 
где х; — точка релевантных данных, которые оцениваются, х(о) — опорная величина 
номинальных данных, с- масштаб естественной вариации, [{ — пороговый параметр. 
Для устойчивой оценки местоположения наблюдаемых данных и их разброса часто 


рекомендуют использовать величину медианы и абсолютное отклонение от медианы 
(МАР) [6]: 


: х и +Х , 
тефап(Х у) = В ГУньм 


5 = МАБСХ у) =1.4826х тей а! х, — тефаи(Х,)| р, 
1=... 


Л=Ь....М№ 
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где [®| - функция округления в меныпую сторону к самому близкому целому 
числу, и Л7.^...., Хмм — ранжирование выборки Х». Поправочный коэффициент 1.4826 
выбран с целью, чтобы 5 совпадало со среднеквадратическим отклонением, когда Ху 
является гауссовым распределением. Основным принципом идентификатора Хампеля 
является замена х; = теФап(Ху) и © = 5. Идентификатор Хампеля эффективен, когда 
данные не изменяются во времени. На практике последовательность данных, например, 
временной мониторинг внутричерепного давления, всегда искажена дрейфом базовой 
линии или другими динамическими событиями. Для того чтобы приспособиться к 
динамическим условиям, идентификатор Хампеля применяется в скользящем окне. 
Это подход известен как фильтр Хампеля. 

Для любого целого числа К скользящее окно данных определяется как 

С 

и данные упорядочиваются по рангам, чтобы получить последовательность 


рава х 


(0) (К) 


2 
где хо) центральный элемент в этом ранжированном списке. Следовательно, 
фильтр Хампеля определяется, как: 


= жа <, 


ИЕ 
х 


в №м= жа >15. 


Как только выброс обнаружен, его величина заменяется медианой в этой точке. Поро- 
говую величину / для обнаружения выброса принимают равной {1 = 3 для устойчивости 
к ошибкам. Отметим, что оба экстремальных предела порогового параметра пред- 
ставляют важные особые случаи: когда 1 = 0, мы возвращаемся к медианному фильтру, 
поскольку у; = хо) всегда поддерживается в этом случае, и в пределе, когда (>00 мы 
приближаемся к фильтру идентичности у; = х;. Фильтрация отсутствует, так как в 
этом случае медианный критерий выбора никогда не может удовлетворяться. 


2 Калмановская фильтрация 


После обнаружения в онлайновом режиме с помощью фильтра Хампеля выбросов 
в наблюдаемом временном ряде данных и замены их медианой данных в скользящем 
окне эти точки сигнала подлежат восстановлению. Важным обстоятельством после 
удаления выбросов является то, что распределение измеряемых данных можно считать 
гауссовым распределением и использовать байесовское вероятностное оценивание 
сигнала из зашумленных наблюдений. В этой методологии задача оценивания рас- 
сматривается в виде последовательной двухэтапной процедуры прогноза и коррекции, 
основанной на предположении аддитивной модели линейных уравнений процессов 
оценивания и наблюдения, дополненных белыми гауссовыми шумами. В этом случае 
распределения вероятностей оценок также оказываются гауссовыми и решение задачи 
сводится к рекуррентной системе уравнений. Практическая реализация этого подхода 
известна как фильтр Калмана, в котором стадия оценивания и прогноза, основанная 
на марковском процессе, является моделью обновления во времени и используется 
для предсказания следующего состояния системы по данным текущего состояния. 
Стадия коррекции использует независимую модель измерения для коррекции пред- 
сказанной оценки по предоставленным данным наблюдения. Фильтр Калмана также 
включает динамический фактор усиления в уравнениях временных обновлений, кото- 
рый автоматически обновляется для точного определения, как много «доверия» вве- 
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дено в наблюдаемые данные. Поэтому, если наблюдения, как считают, очень зашум- 
лены, больше «доверия» размещено в прогнозирующую стадию, и таким образом 
стадия коррекции оказывает небольшое влияние на сигнал. Для реализации фильтра 
Калмана необходима также начальная оценка состояния системы (значения переменных 
состояния) и ковариации ошибки оценки. Уравнения фильтра Калмана детализированы 
в книге [7]. 

Общим подходом для распознавания и прогноза сигнала является модель прос- 
транства состояний, которая может использоваться для описания большого количества 
различных моделей. Одной из таких моделей является адаптивная авторегрессионная 
ААК(р) модель р-го порядка. Параметры ААК могут использоваться для описания 
изменяющихся во времени данных, охватывая второй порядок статистических моментов. 
В этой модели не требуется почти никакого априорного знания, порядок модели р не 
очень важен и, так как он является одиночным коэффициентом, он может быть легко 
оптимизирован. Кроме того, нет необходимости в выборе особо дорогого алгоритма. 
Параметры ААК обеспечивают простой и робастный подход и, следовательно, хорошую 
отправную точку для адаптивного выявления характерных особенностей. Одномерная 
авторегрессионная (АК) модель, описывается следующим уравнением: 


р 
х(п) = У ‘а, х(п-Ю)+е(п), 
= 
где х(и)- анализируемый физиологический сигнал (например, сигнал ВЧД) в 


момент времени и, {а ра — параметры модели, НС ЮР} — задержанные выборки 


сигнала с процессом обновления е„ = и, = 0, о) с нулевым средним значением и дис- 
персией о”. Чтобы выполнить оценку ААВ параметров фильтром Калмана, модель ААК 
должна подходящим способом быть приспособлена к модели пространства состояний. 
С этой целью параметры АК(р) становятся векторами состояния &4 = ах = [а14 ``, ар] я 
Решение проблем выбора начальных значений, коэффициентов обновления и порядка 
модели фильтра Калмана для оценки ААК параметров, а также его применения для 
обработки сигналов электроэнцефалограмм (ЭЭГ) можно найти в работе [8]. Ранее 
мы использовали алгоритмы калмановской фильтрации и калмановского сглаживания 
для адаптивной оценки параметров авторегрессионной модели при изучении морфо- 
логии волн ВЧД при острой травме головного мозга [9]. 

Применение авторегрессионной модели, в частности, авторегрессионой модели 
проинтегрированного скользящего среднего (АЕВЛМА) с выбором порядка на основе 
автокорреляционной (АСЕ) и частной автокорреляционной (РАСР) функций для пред- 
сказывания непрерывных трендов сигнала ВЧД было предложено в работе [10]. Было 
показано, что точность предсказания ВЧД значительно лучше, кода порядок модели 
АЕПМА оценивается АСЕ и РАСЕ, чем с помощью информационного критерия Акаике. 
Этими же авторами был предложен алгоритм, использующий авторегрессионную 
модель, для восстановления сигналов ВЧД после удаления артефактов и выбросов с 
помощью эмпирического метода декомпозиции [11-12]. 


Результаты и обсуждение 


Инвазивные измерения внутричерепного давления (ВЧД) осуществлялись с помо- 
щью универсального сетевого мониторинго-диагностического комплекса для управ- 
ляемой терапии больных. Составляющими звеньями этого комплекса были: прибор 
мониторинга ВЧД, основанный на жидкостно-заполненной системе катетер-тензометр 
с измерительным модулем М3500Р (УкрЦЕНДИСИ НАН Украины, г. Киев) [13], и 
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Со4дтап ГСР Ехргез$ монитор с тензометрическими микродатчиками давления [14]. 
Системное программное обеспечение позволяет адресное обращение к измерительным 
модулям, введение информации к компьютеру, накопление и математическую обра- 
ботку данных и их графическое изображение в течение мониторинга несколько суток 
непрерывно. Мониторинг ВЧД проводился у 117 пострадавших в возрасте от 18 до 
62 лет с тяжелой черепно-мозговой травмой (глубина коматозного состояния по 
шкале ком Глазго менее 8 баллов), находящихся на лечении в отделении нейрореа- 
нимации Донецкого областного клинического территориального медицинского объе- 
динения (ДОКТМО). 

На рис. 1 а показан типичный пример сигнала ВЧД, неинвазивно регистриру- 
емого в клинических условиях нейрореанимации, который искажен множественными 
выбросами из-за неконтролируемого движения пациентом головой. Процесс неавто- 
номной предварительной обработки сигнала фильтром Хампеля показан на рис. 2, 
где вместе с измеренным ВЧД сигналом приведены результаты обработки медианным 
фильтром, верхняя и нижняя пороговые величины обнаружения выбросов, обнару- 
женные выбросы и конечный результат фильтрации. К обрабатываемым данным сиг- 
нала применено скользящее окно, сформированное последовательностью данных {хи} 
длинной К точек. Чувствительность фильтра Хампеля контролируется двумя пара- 
метрами — порогом обнаружения выбросов — [и К — шириной окна. Как упоминалось 
выше, пороговая точка средней статистической величины или отношение числа выб- 
росов к общему количеству выборок в последовательности данных, для которой средняя 
статистическая величина становится выбросом, составляет 50%. Это обстоятельство 
предполагает, что ширина окна К должна быть, по крайней мере, вдвое шире отдель- 
ного выброса. Если в пределах ширины окна присутствуют два или больше близких 
выброса фильтр обнаружит их присутствие, если их суммарная ширина вдвое меньше 
ширины окна. 


Опота! здпа! 
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Рисунок 1 — Сравнение сигналов 138-секундного фрагмента мониторинга ВЧД 
сигнала у пациента Кос-ко: а) до и Б) после обработки фильтром Хампеля в 
неавтономном режиме. Точки выбросов заменены значением медианы в этих точках 
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В работе [15] исследовалась возможность применения фильтра Хампеля для 
удаления выбросов в сигналах электроэнцефалограмм, где было найдено, что К`= 80 
измерительных точек и величина порога -5 — 6 были подходящими для всех субъектов. 

Для обнаружения выбросов в сигнале ВЧД, показанных на рис. | и рис. 2, мы 
использовали ширину окна К = 210 точек и пороговую величину равной / = 3. Обра- 
ботанный фильтром Хампеля сигнал ВЧД показан на рис. 1Ъ, где обнаруженные точки 
выбросов заменены на значение медианы в этих точках. 

Задача второго (заключительного) этапа обработки состояла в фильтровании 
полученного после удаления выбросов сигнала — медианные значения следует заменить 
на «хорошие» данные, которые с наибольшей вероятностью будут соответствовать 
структуре исследуемого сигнала. Для этой цели использовалась ААЁК модель, пара- 
метры модели определялись с помощью алгоритма Калмана и последующим КТ$ 
сглаживанием. Фильтр Калмана представляет собой универсальный инструмент для 
онлайновой оценки коэффициентов модели. Однако производительность фильтра 
Калмана существенно зависит от его инициализации, выбора коэффициентов обнов- 
ления и порядка модели, которые обычно неизвестны. Конечно, их можно принять 
навскидку, например, нулевыми, но такой подход приводит к длительному начальному 
переходному эффекту. Кроме того, нулевая установка начальных значений парамет- 
ров ААК может иметь некоторое вредное влияние на результат. Основываясь на неко- 
торых предварительных экспериментах, в работе [8] было показано, что в модели 
пространства состояний для запуска алгоритма рекомендуется ковариационная матрица 
начального состояния в виде диагональной матрицы, взвешенной дисперсией преды- 
дущих данных и умноженной на коэффициент д, который может принимать либо 
очень малые значения, либо очень большие О) = до. В нашем случае мы приняли д= 10°. 
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Рисунок 2 — Предварительная обработка фильтром Хампеля фрагмента сигнала 
мониторинга ВЧД, полученного в клинике у пациента Кос-ко. 
Ширина скользящего окна взята размером К`= 120 точек. Штрихпунктирная линия 
показывает параметры прогноза опорной величины номинальных данных 
хо) = шеФап(х). Пунктирные линии показывают прогнозируемый порог выбросов 
(хо) + 35). Точки обнаруженных выбросов показаны квадратом (5). 
Точечной линией показан сигнал ВЧД, очищенный от выбросов 
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Выбор коэффициента обновления подробно обсуждался в работе [16] и было 
установлено, что в случае ААК оценки сигналов ЭЭГ оптимальный выбор коэффициен- 
та обновления является компромиссным решением между скоростью адаптации и то- 
чностью оценки. С этой целью, ковариационная матрица шума процесса И» и дисперсия 
процесса обновления Их в зависимости от коэффициента обновления ИС были взяты в 
виде: 

И’, = [-ИСлтасе( Он, 
РЕ ЬС, 

где коэффициент обновления в нашей работе мы приняли равным ОС = 48 [16]. 

Оптимальный порядок модели мы определяли, используя информационный 
критерий Акаике (АПС) и байесовский информационный критерий (ВС). Результаты та- 
кой оценки представлены на рис.3Ъ. К сожалению, порядок модели, полученный этими 
подходами, не был обязательно лучшим для обработки сигналов ВЧД других пациентов 
или разных фрагментов длительного мониторинга одного пациента. Наш опыт показы- 
вает, для устойчивой работы алгоритма ВТЗ сглаживания часто требуются намного 
меньшие порядки, чем определяют критерии АГС и В]С. В основном, мы использовали 
порядок 5 или 6. Из доступной нам литературы, мы полагаем, что вопрос оптимальных 
начальных параметров ААК оценки остается все еще открытым. 
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Рисунок 3 — Результат онлайновой двухэтапной обработки ВЧД сигнала, 
загрязненного множеством выбросов — артефактов движения, возникших 
при мониторинге пациента К-ко в клинических условиях нейрореанимации: 
(а) загрязненный ВЧД сигнал; (Б) очищенный сигнал после обнаружения и замены 
выбросов фильтром Хампеля, прогноза замененных данных фильтром Калмана 
и окончательным КТЗ сглаживанием 


Для этого случая наилучшая оценка р = 13 порядка ААК модели была опреде- 
лена ВГС критерием. 
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Выводы 


В этой статье представлена методика неавтономного фильтрования физиологи- 
ческих сигналов. Исследование сосредоточено на обнаружении выбросов и очистке 
сигналов ВЧД мониторинга. Удаление выбросов является необходимым шагом пред- 
варительной обработки жизненно важных медицинских сигналов, чтобы усилить их 
характеристики для дальнейшего диагноза. Данная схема фильтрования осуществля- 
ется в два этапа — сначала используется робастный фильтр Хампеля со скользящим 
окном для обнаружения и удаления выбросов, затем калмановская фильтрация и кал- 
мановское сглаживание для оценки и восстановления ВЧД сигнала. Эффективность 
предложенного метода удаления артефакта экспериментально оправдана, основываясь 
на сигналах ВЧД мониторинга 117 больных в возрасте от 18 до 62 лет с тяжелой черепно- 
мозговой травмой; выбросы и шум удаляются без искажений оригинального ВЧД 
сигнала. Этот метод может также использоваться для фильтрации любых нестацио- 
нарных техногенных сигналов. 
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ОпПпе оиШегу айесйоп апа с(еатие 


орттгастата ргеззиге топпоттв яепа[5 

Те уа5Ё та]отЦу ог тезеагспез ш фе 41юИа| даа ргосеззше геЁег ю Ше 512па1$ оЁ 
апгоровешс опеш, ш \Лмеь ргоегл$ оЁ по15е за БИщу апа по15е НЦегие оЁ по15у 51е1а[$ 
сепегае4 Бу Четепе теап$ оР сотититсайоп ап4 БгоадсазНпо аге зо]уе4. АЕ ]а5{ 4есадез, 
езреслаПу абгоа4, Вауе Бееп ицепзтуе[у 4еуе!ореЯ сотршег-Базе4 ргосеззше оЁ рНуз10- 
1оэ1са| $1епа15 ш теса!| аррИсаНоп$. ТВе рВуз1о]оглса1 $1епа15 апа]уз15 Уаз плисВ Ваг4дег 
Фап ап@горосетс 51юпа[5. ТВаф 15зие дие Бу Ше мо Рсюгз: рнуз10]ое1са| ргосеззез Ваз 
$12па[$ ул 1о\/ атрПаде, ап4 езе $1епа15 аге шНегепИу поп-%айопагу ап4 поп-Ппеаг. 
п ада оп, Фезе $1опа[5 аге соматштаеа ВеауПу зюсвазИс сотропеп, Гог ехатр/е, Пе 
$1епа| топНоние оЁ шёасгаша| ргеззиге (ТСР), гезИпе ш Меиго Пуцепууе Саге ЧпИз 
оНеп соша1$ а 1агое аточп: оЁ по15е ап оиЙетгз. ТБезе агИас{$ по ошу ЯтесИу |еа4 0 
Ёа]5е а[агтз ш ашютайс а|е{ зузетз оЁГ топйогпе-ФаепозИс сотр!ех Юг соптоПед 
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{егару раНепь, ап Шеу а[50 зеуегейу соматштае Фе сБагащен$Нсз оР Фе ипдейуте 
31епа|, уЛасВ такКез ассигайе ГогесазИпе о{ зесопдагу Бгат Чатасез саизе4 Бу шёгасгата! 
Нуремепз1оп ппро$$1Ые. 

Аррагепйу ех15{5, {гее Баз1с НЦег Рог даесНоп ап4 гетоуа! оР оч ег$ Нот ргосез$ 
ава зепез, УЛмсВ Ва тодегитаноп ш гесеп{ \уогк$. ТЫ$ 15 фе тефап НШег, с1озейу те 
{0 И ап ао аз зпаре - Натре! НЦег ап Ве НЦег-Феапег, 4ез1юптед Бу Маги апа 
Тротзоп (МТ НИег-с‹еапег). Ш Ве Вгз{ хо НШег$ сап \’отК ш Фе оп-Ппе тоде, ве МР 
ВНЦег-Сеапег — уогК$ Бу еззепнаПу аиюпотоп$ зсВете, умей 15 Базе4 оп а рге-езитае4 
то оЁ Фе Ите зепез. ТВезе йгее Баз1с НШегз ууеге {е5е4 ог Фе 14епиЙсаНоп оЁ гее 
Чегет Чупапис зузеплз, опе оЁ \ЙсВ уаз нежей у геа! оп Йег по1зу з1юпа[ оЁ Фе 
Бейсорег Ве1о{ ап@ апз]е соплтап4$ ап@ Фе о ег {о тоде| зузетз — Ппеаг то4е]| ап4 
попПпеаг зесоп4-ог4ег Уо[ета то4е| у а9Фуе Саяз$ап по1зе ап оиШег$, рас 
аге Кпо\л Юг ассига{е гези5. [шп Гасё, Феге 15 по еЙесНуе итуегза| с1еаппе НШ ег даа. 

ТЫ$ ага‹е оНегз ап еРйсеп{ оп-Ппе \о-“ер риййсаНоп тефо4 оЁ рБуз10]о51са1 
$1юпа[$ Базед оп Натре| 14епиНег апа Ка|тап НПегто. шшаПу, сПиса! теазигетепЕ оЁ 
ТСР япа| ипдегоо ргейшеантеп ш УЛисВ Фе 14епййсайоп оЁ 51юпа| даа збтасвге аге геаП7е4, 
ап фе езитаеа [еуе] оЁ по15е оч Иег$ аге гетоуеа Бу Фе Натре| НИег гоба$. Ош@ег 
ро аге гер|асед Бу фе те ап уаше а{ Фозе рош6. Мех, Фе согесНоп оЁ Фезе рошк, 
п015е гетотуа| ап еуа[ааНоп оЁ 1е $1епа| Базе оп адарнуе амюгеотезяуе (ААВ) то4е! 
изше Фе Ка|тап НШег ап4 5 аззостаеа ВТЗ (Каисв-Типе-З6леБе!) зтоо ше НЦег аге 
регогтеа. Рторозе4 Нет» тео сизюпита е рагатеегз аге Фе Ва№-ул4 тоуше 
улпдоу,, фе Шгезро!4 Гог даесипе оиегз ап4 ог4ег ААК тоде|. Те еНесйуепез$ оЁ 
{Ве ргорозе4 тео 15 а5ИНеа Бу Ше ехрегипеща! агас( гетоуа1, Базе оп Ше $1епа15 
оЁ [СР топйопие аё 117 райет асе 18 ю 62 уеагз УИ 5еуеге гаитайс Бгаш шуагу, 
ош ег ап по15е аге гетоуе улой @Язюгйоп оЁ Фе опрта| СР з1епа[. ТЬ1$ шефо4 
сап а15о Бе изе4 1ю НЦег оц апу поп-$айопагу ап@горозешс з12па|5. 


Статья поступила в редакцию 26.04.2013. 
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